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Abstract 



^ ■ The analysis of shape transitions in Nd isotopes, based on the framework of relativistic energy 

OO . density functionals and restricted to axially symmetric shapes in Ref. [22], is extended to the re- 

gion Z = 60, 62, 64 with N ~ 90, and includes both (3 and 7 deformations. Collective excitation 
■ spectra and transition probabilities are calculated starting from a five-dimensional Hamiltonian 

. for quadrupole vibrational and rotational degrees of freedom, with parameters determined by con- 

. ^ strained self-consistent relativistic mean- field calculations for triaxial shapes. The results reproduce 

X 

^ ' available data, and show that there is an abrupt change of structure at = 90 that can be approx- 

^ '. 

imately characterized by the X(5) analytic solution at the critical point of the first-order quantum 
phase transition between spherical and axially deformed shapes. 
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I. INTRODUCTION 



The evolution of shell structures governs the variation of ground-state nuclear shapes 
along isotopic and isotonic chains. Nuclear structure explores a variety of phenomena related 
to structural evolution including, for instance, the reduction of spherical shell gaps and 
modifications of magic numbers in nuclei far from stability, occurrence of islands of inversion 
and coexistence of shapes with different deformations. As the number of nuclcons changes 
from nucleus to nucleus, in general one observes a gradual evolution of different shapes 
- spherical, axially deformed, shapes that are soft with respect to triaxial deformations. 
These shape transitions reflect the underlying modifications of single-nucleon shell structure 
and interactions between valence nucleons. An especially interesting feature is the possible 
occurrence of shape phase transitions and critical-point phenomena for particular values of 
the number of protons and neutrons. Phase transitions in the equilibrium shapes of nuclei 
correspond to first- and second-order quantum phase transitions (QPT) between competing 
ground-state phases induced by variation of a non-thermal control parameter (number of 
nucleons) at zero temperature [1-3]. 

Nuclear quantum phase transitions have been the subject of extensive experimental and 
theoretical studies during the last decade. For recent reviews and an exhaustive bibliography 
we refer the reader to Refs. [4-8]. Even though phase transitions in finite systems can 
only be defined in the classical limit in which the number of constituents tends to infinity, 
i.e. the transition is actually smoothed out in finite systems, there are nevertheless clear 
experimental signatures of abrupt changes in structure properties of finite nuclei with the 
addition or subtraction of only few nucleons. This is another distinct feature of QPT in 
atomic nuclei, i.e. the physical control parameter - number of nucleons, can only take 
discrete integer values. Expectation values of suitably chosen operators, that as observables 
characterize the state of a nuclear system, can be used as order parameters [9]. 

Theoretical studies of nuclear QPT are typically based on macroscopic geometric mod- 
els of nuclear shapes and potentials [8] and/or semi- microscopic algebraic models [7]. In 
the geometric framework QPTs are analyzed in terms of a Bohr collective Hamiltonian 
for shape variables, and can be related to the concept of critical-point symmetries that 
provide parameter-independent predictions for excitation spectra and electromagnetic tran- 
sition rates for nuclei at the phase transition point. Alternatively, in the algebraic approach 
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(e.g. the interacting boson model (IBM)) different shapes coincide with particular dynamical 
symmetries of some algebraic structure, and the phase transition occurs when these symme- 
tries are broken in a specific way. The two frameworks can be related, for instance, by using 
the coherent state formalism [10, 11] which allows to establish a correspondence between the 
symmetry limits of the algebraic Hamiltonian and energy functionals expressed in terms of 
collective shape variables. In both approaches, geometric and algebraic, the description of 
QPT is based on model specific Hamiltonians which by construction describe shape changes. 
A shape phase transition is then accessed by variation of a control parameter. 

The two best studied classes of nuclear shape phase transitions, both theoretically and 
experimentally, correspond to a second-order QPT between spherical and 7-soft shapes 
[12, 13], and a first-order QPT between spherical and axially deformed shapes [14, 15]. The 
former is a phase transition in one degree of freedom - the axial deformation (3 and, in the 
IBM language, represents a transition between the U{h) and 0(6) dynamical symmetries in 
the limit of large boson number. The critical point of phase transition can also be related to a 
dynamical symmetry: in this case E(5) [12], and the experimental realization of this critical- 
point symmetry was first identified in the spectrum of ^^"^Ba [13]. The second type of shape 
transitions, between spherical and axially deformed nuclei, is more commonly encountered 
and involves two degrees of freedom - the deformations j3 and 7. The critical point of this 
phase transition, denoted X(5), does not correspond to a dynamical symmetry in the usual 
sense [14]. Nevertheless, for the particular case in which the (3 and 7 degrees of freedom 
are decoupled, an approximate analytic solution at the critical point of phase transition can 
be expressed in terms of zeros of Bessel functions of irrational order [14] . Evidence for the 
empirical realization of X(5) phase transition was first reported for ^^^Sm and other iV = 90 
isotones in Rcf. [15]. 

Even though phenomenological approaches to nuclear QPT have been very successful, 
and the predicted isotopic trends for various observables are in very good agreement with 
data, it would clearly be desirable to have a fully microscopic description of shape phase 
transitions, starting from nucleonic degrees of freedom. This is especially important in 

view of the fact that the physical control parameter in nuclear QPT is the actual number 
of miclcons in a nucleus, rather than a strength parameter of a model specific, Ising-type 
Hamiltonian. 

Several recent studies [16-21] have reported microscopic calculations of potential energy 
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surfaces as functions of quadrupole deformation, for a number of isotopic chains in which the 
occurrence of shape phase transitions had been predicted. Potentials calculated with a con- 
straint on the mass quadrupole moment display shape transitions from spherical to deformed 
configurations. In Ref. [22] we analyzed shape transitions in Nd isotopes, using a microscopic 
model based on relativistic energy density functionals (EDF). Starting from constrained 
self-consistent mean-field calculations of potential energy curves, the generator coordinate 
method (GCM) [23] was used to perform configuration mixing of angular-momentum and 
particle-number projected relativistic wave functions. It was shown that the microscopic 
framework based on universal EDFs, adjusted to nuclear ground-state properties, and ex- 
tended to take into account correlations related to symmetry restoration and fluctuations 
of collective variables, describes not only general features of shape transitions, but also the 
unique behavior of excitation spectra and transition rates at the critical point of quantum 
shape phase transition. In particular, the particle-number projected GCM spectra, intra- 
band and interband B(E2) values for I48,i50,i52]yj^ were found in excellent agreement with 
data, and close to the X(5)-model predictions for ^^'^Nd. The self-consistent GCM calculation 
based on the relativistic density functional PC-Fl [24], predicted the shape phase transition 
precisely at the isotope ^^°Nd, in agreement with empirical evidence for the realization of 
X(5) in the N — 90 rare-earth isotones [5, 6]. 

X(5) denotes a particular model for a first-order QPT between spherical and axially 
deformed shapes, based on the assumption of a separable potential V{P, 7) = Vp{P) + Vy{'y). 
By neglecting the small barrier between the competing spherical and deformed minima, i.e. 
taking an infinite square well in the variable (3, and assuming an axially stabilized potential 
V^{'j), an approximate analytic solution is obtained at the point of phase transition [14]. The 
microscopic axially symmetric potential for ^^°Nd, calculated with self-consistent mean-field 
models (cf. Fig. 1 in Ref. [22]) is, of course, more realistic than the infinite square well 
considered by lachello in Ref. [14]. First, although it displays a wide fiat minimum in /3 on 
the prolate side, the fiat bottom of the potential does not start at /3 = 0, i.e. the coexisting 
shapes do not include the spherical configuration. In this respect the microscopic potential 
is closer to the generalization of the axially symmetric X(5) solution for the collective Bohr 
Hamiltonian to the transition path between X(5) and the rigid-rotor limit, represented by 
infinite square-well potentials over a confined range of values (3m > /?m ^ (confined 
soft potential, and /3m denote the positions of the inner and outer walls on the /3 axis. 
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respectively) [25]. Second, the microscopic potential well has softer walls compared to the 
X(5) square well. The unrealistic features introduced by the hard wall of the square well 
were analyzed in Ref. [26], and it was noted that the compression of the wave function 
against the well wall induces effects in spectra approximating rigid (3 deformation. 

More importantly, the GCM calculation of Ref. [22] was restricted to axially symmetric 
shapes (7 = 0). We note that an exactly separable 7-rigid version (with 7 = 0) of the 
X(5) model has been constructed in Ref. [27]. Although the original X(5) solution relied 
on an approximate separation of the potential in the variables (5 and 7 [14], the exact 
diagonalization of the Bohr Hamiltonian with a /3-soft, axially stabilized potential, carried 
out in Ref. [26], has shown that many properties of the solution are dominated by /? — 7 
coupling induced by the kinetic energy operator. Band-head excitation energies, energy 
spacings within the bands, and transition strengths are strongly dependent on the 7 stiffness 
of the potential. The importance of the exphcit treatment of the triaxial degree of freedom, 
i.e. inclusion of /3 — 7 coupling, was also emphasized in two recent studies [20, 21] that used 
the self-consistent Hartree-Fock-Bogoliubov model, based on the finite-range and density- 
dependent Gogny interaction, to generate potential energy surfaces in the /? — 7 plane. 
In addition, in Ref. [20] excitation spectra of Nd isotopes were calculated using the GCM 
with particle-number and angular-momentum projected axially symmetric wave functions. 
However, while GCM configuration mixing of axially symmetric mean-field states is routinely 
performed in studies of collective excitation spectra, the application of this method to triaxial 
shapes is a much more difficult problem. Only very recently a model has been introduced [28] , 
based on triaxial mean-field states, projected on particle number and angular momentum and 
mixed by the generator coordinate method. The numerical implementation of the model is 
rather complex, and applications to medium-heavy and heavy nuclei are still computationally 
too demanding. An alternative approach is the explicit solution of the collective Hamiltonian 
in five dimensions, with deformation-dependent parameters determined from microscopic 
self-consistent mean-field calculations. This approach will be used in the present study. 

In this work we therefore extend the analysis of shape transitions in Nd isotopes of 
Ref. [22], to the region Z = 60, 62, 64 with ^ 90, and include the explicit treatment of both 
(3 and 7 degrees of freedom. Collective excitation spectra and transition probabilities are 
calculated starting from a five-dimensional collective Hamiltonian for quadrupole vibrational 
and rotational degrees of freedom, with parameters determined by constrained self-consistent 
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relativistic mean-field calculations for triaxial shapes. 

An outline of the theoretical framework is included in Sec. II: the method of solution of 
the eigenvalue problem of the general collective Hamiltonian, and the calculation of the mass 
parameters, moments of inertia, and zero-point energy corrections. In Sec. Ill we present a 
detailed study of shape transitions in the N ^ 90 rare-earth nuclei. Theoretical results are 
compared with experimental excitation spectra and transition rates, as well as with the X(5) 
approximate solution at the point of first-order quantum phase transition. Sec. IV presents 
a summary and an outlook for future studies. 



II. COLLECTIVE HAMILTONIAN IN FIVE DIMENSIONS 



In Ref. [29] we have developed a model for the solution of the eigenvalue problem of a 
five-dimensional collective Hamiltonian for quadrupole vibrational and rotational degrees of 
freedom, with parameters determined by constrained self-consistent relativistic mean-field 
calculations for triaxial shapes. The five quadrupole collective coordinates are parameterized 
in terms of two deformation parameters (3 and 7, and three Euler angles (0, 9, ip) = ft, 
which define the orientation of the intrinsic principal axes in the laboratory frame. The 
collective Hamiltonian can be written in the form: 



H — Tvib + Trot + v. 



'coll 



(1) 



with the vibrational kinetic energy: 



vib 



+ 



1 



w 07 



/3 sin 37 
and rotational kinetic energy: 
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(2) 



(3) 



V^ou is the collective potential. denotes the components of the angular momentum in 
the body-fixed frame of a nucleus, and the mass parameters Bpp, Bpj, B^^, as well as the 
moments of inertia Ik, depend on the quadrupole deformation variables (3 and 7: 



Ik = 4Sfe/32sin2(7 - 2A;7r/3) 



(4) 
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Two additional quantities that appear in the expression for the vibrational energy: r = 
S1S2-B3, and w — BppB^j — determine the volume element in the collective space. The 
Hamiltonian Eq. (1) describes quadrupole vibrations, rotations, and the coupling of these 
collective modes. The corresponding eigenvalue problem is solved using an expansion of 
eigenf unctions in terms of a complete set of basis functions that depend on the deformation 
variables /? and 7, and the Euler angles 0, 9 and ip [30]. 

The dynamics of the collective Hamiltonian is governed by the seven functions of the 
intrinsic deformations (3 and 7: the collective potential, the three mass parameters: Bpis, 
Bi3^, B^^, and the three moments of inertia Ik- These functions are determined by the 
choice of a particular microscopic nuclear energy density functional or effective interaction. 
As in our previous study of axial shape transitions in Nd isotopes [22], also in this work we 
use the relativistic functional PC-Fl (point-coupling Lagrangian) [24] in the particle- hole 
channel, and a density-independent 5-force is the effective interaction in the particle-particle 
channel. The parameters of the PC-Fl functional and the pairing strength constants T4 and 
Vp have been adjusted simultaneously to ground-state observables (binding energies, charge 
and diffraction radii, surface thickness and pairing gaps) of spherical nuclei [24], with pairing 
correlations treated in the BCS approximation. 

The map of the energy surface as function of the quadrupole deformation is obtained by 
imposing constraints on the axial and triaxial mass quadrupole moments. The method of 
quadratic constraints uses an unrestricted variation of the function 

{H) + J2 {{Q2,) - q2,y , (5) 

M=0,2 

where (H) is the total energy, and (02^) denotes the expectation value of the mass 
quadrupole operator: 

(520 = - - and Q22 = - (6) 

is the constrained value of the multipole moment, and C21X the corresponding stiffness 
constant [23]. 

The moments of inertia are calculated according to the Inglis-Belyaev formula: [31, 32] 

^=1,2,3, (7) 
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where k denotes the axis of rotation, and the summation runs over the proton and neu- 
tron quasiparticle states. The quasiparticle energies Ei, occupation probabihties Vi, and 
single-nucleon wave functions V'i are determined by solutions of the constrained relativis- 

tic mean-field (RMF) plus BCS equations. The mass parameters associated with the two 
quadrupole collective coordinates go = (Q20) and ^2 = {Q22) are also calculated in the 
cranking approximation [33] 

(8) 
with 

M(n),i.u{qQ, q2) = 2^ (E- + EY ^'^^ ' 

i,j v * ]/ 

The collective energy surface includes the energy of zero-point motion (ZPE), that corre- 
sponds to a superposition of zero-point motion of individual nucleons in the single-nucleon 
potential. The ZPE correction depends on the deformation, and includes terms originating 
from the vibrational and rotational kinetic energy, and a contribution of potential energy 

^V{qo, q2) = AKib(go, q2) + AKot(go, ^2) + AVpot(go, ^2) ■ (10) 

The latter is much smaller than the kinetic energy contribution, and is usually neglected. 
The vibrational and rotational ZPE are calculated in the cranking approximation [33], i.e. 
on the same level of approximation as the mass parameters and the moments of inertia. 
The potential Kou in the collective Hamiltonian Eq. (1) is obtained by subtracting the ZPE 
corrections from the total energy that corresponds to the solution of constrained RMF-I-BCS 
equations (cf. Eq (5) ), at each point on the triaxial deformation plane: 

Koii(go, 52) = -Etot(?o, 52) - AKib(?o, 52) - AVrotiqo, 52) ■ (11) 

III. SHAPE PHASE TRANSITIONS IN THE N a 90 REGION 

The present study considers the Nd, Sm and Gd nuclei with N ^ 90. These isotopic 
chains display transitions from spherical to deformed shapes, and experimental evidence has 
recently been reported [15, 34, 35] for the N ^ 90 isotones: ^^°Nd, ^^^Sm, ^^"^Gd lying close 
to the point of first-order QPT. We will first extend our analysis of Ref. [22] of the Nd 
isotopic chain to explicitly include the triaxial degree of freedom. The microscopic picture 



Bf,u{qo:q2) = Y 



•^(i;-^(3)-M(i; 
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of shape phase transitions at = 90 will be analyzed in detail for the Nd nuclei, and in the 
remainder of this section the most interesting results for the Sm and Gd isotopic chains will 
be presented. 

In Fig. 1 we display the self-consistent RMF+BCS triaxial quadrupole binding energy 
maps of the even-even ^'^^~^^"^Nd in the /3 — 7 plane (0 < 7 < 60*^), obtained by imposing 
constraints on expectation values of the quadrupole moments {Q20) and {Q22) (cf. Eq. (5)). 
Filled circle symbols denote absolute minima; all energies are normalized with respect to 
the binding energy of the absolute minimum. The contours join points on the surface with 
the same energy. The energy maps nicely illustrate the gradual increase of deformation of 
the prolate minimum with increasing number of neutrons, from the spherical ^^"^Nd to the 
strongly deformed ^^"^Nd, and the evolution of the 7-depcndcncc of the potentials. With 
increasing neutron number the Nd isotopes remain prolate deformed in the lowest state, 
i.e. the shape evolution corresponds, in the language of the interacting boson model, to a 
transition between the U(5) and SU(3) limits of the Casten symmetry triangle [5, 6]. One 
notes, however, the appearance of oblate minima in ^^^Nd and ^^^Nd. An important feature 
of the energy maps shown in Fig. 1 is the extended flat minimum around j3 ~ 0.3 in ^^°Nd. 
This is seen even more clearly in the axial maps of Fig. 2, where we plot the binding energies 
as functions of the axial deformation parameter (3. Energies are normalized with respect 
to the binding energy of the absolute minimum for a given isotope. Negative values of j3 
correspond to the oblate (/3 > 0, 7 = 60°) axis on the /3 — 7 plane. The prolate minimum 
gradually shifts to larger deformation and saturates at /3 0.35 in ^^^Nd and ^^^Nd. It 
is interesting that for these two nuclei the calculation also predicts the same deformation 
j3 ~ —0.3 of the oblate minima. The difference is in the barriers separating the prolate 
and oblate minima: 6 MeV for ^^^Nd, and more than 7 MeV for ^^"^Nd. For ^^''Nd, on the 
other hand, the flat prolate minimum extends in the interval 0.2 < f3 < 0.4, whereas the 
oblate minimum at /3 ~ —0.2 is not a true minimum, but rather a saddle point in the — 7 
plane (cf. Fig. 1). This is illustrated in Fig. 3, where we plot the binding energy curves of 
^^•^Nd as functions of the deformation parameter 7, for two values of the axial deformation 
j3 = 0.2 and 0.25. In the region of the flat prolate minimum the potential displays a parabolic 
dependence on 7 for 7 < 30°, but flattens with increasing 7 because, of course, it must be 
periodic with a period 27r/3. The difference between the two curves plotted in Fig. 3 shows 
that the 7-potential is not completely independent of /3, as it was originally assumed for 
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the X(5) model in Ref. [14]. We note that both in the test of the X(5) for the 7 degree of 
freedom of Ref. [40], as well as in the study of the validity of the approximate separation 
of variables introduced with the X(5) model [26], the potential in 7 was chosen to be a 
harmonic oscillator potential. 

Potential energy surfaces similar to those shown in Fig. 1 were also obtained in two recent 
studies [20, 21] that used the self-consistent Hartree-Fock-Bogoliubov (HFB) model, based 
on the finite-range and density-dependent Gogny interaction. However, compared to the 
results shown in Figs. 1 and 2, the fiat prolate minimum in ^^°Nd is less pronounced in the 
HFB calculation with the Gogny force. By performing the corresponding GCM calculation 
with particle-number and angular-momentum projected axially symmetric wave functions, 
the authors of Ref. [20] show that, rather than ^^°Nd, the excitation energies and transition 
probabilities for ^^*Nd are actually closer to the X(5) model predictions. In fact, they notice 
that intrinsic potentials rather different from the one considered in the original X(5) model 
[14] (wide fiat minimum in (5 starting at /3 = 0), lead to excitation spectra and transition 
rates similar to the X(5) predictions, and thus question the use of the axial deformation (5 
as order parameter of a shape phase transition. 

Starting from constrained self-consistent solutions, i.e. using single-particle wave func- 
tions, occupation probabilities, and quasiparticle energies that correspond to each point on 
the energy surfaces shown in Figs. 1, the parameters that determine the collective Hamilto- 
nian: mass parameters Bpp, Bp^, B^^, three moments of inertia 2^, as well as the zero-point 
energy corrections, are calculated as functions of the deformations P and 7. The diagonal- 
ization of the resulting Hamiltonian yields the excitation energies and the collective wave 
functions: 




(12) 



KeAi 



The angular part corresponds to linear combinations of the Wigner functions 



21 + 1 



(13) 



and the summation in Eq. (12) is over the allowed set of the K values: 




(14) 
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In addition to the yrast ground-state band, in deformed and transitional nuclei excited states 
are usually also assigned to (quasi) (3 and 7 bands. This is done according to the distribution 
of the angular momentum projection K quantum number. Excited states with predominant 
K = 2 components in the wave function are assigned to the 7-band, whereas the /3-band 
comprises states above the yrast characterized by dominant K = Q components. The mixing 
of different intrinsic configurations in the state | al) can be determined from the distribution 
of the projection K of the angular momentum I on the z axis in the body-fixed frame: 

/'7r/3 /»oo 

Nk^^ / 7)r/?1sin37|ci/3o?7, (15) 

Jo Jo 

where the components "ipixiPyl) defined in Eq. (12). When i^' is a good quantum 
number, only one of the integrals Eq. (15) will give a value close to 1. A broader distribution 
of Nk values in the state \al) provides a measure of mixing of intrinsic configurations. 

Before comparing the calculated excitation spectrum of ^^°Nd with data and X(5) model 
predictions, we will consider some characteristic signatures of shape transitions in Nd iso- 
topes. One of the distinct features of shape transitions is the ratio between the excitation 
energies of the first 4"*" and 2"*" states along an isotopic chain. For a transition between spheri- 
cal and axially deformed shapes, this quantity varies from the value Ri/2 = E{4:f) / E{2~1) — 2 
characteristic for a spherical vibrator (U(5) symmetry limit), to R4/2 — 3.33 for a well de- 
formed axial rotor (SU(3) symmetry hmit). In the left panel of Fig. 4 we plot the theoretical 
values of R4/2 for the six Nd isotopes, calculated from the spectrum of the collective Hamil- 
tonian with mass parameters, moments of inertia and potentials determined by the PC-Fl 
density functional, in comparison with experimental values. The calculation reproduces 
in detail the rapid increase of R4/2 from the spherical value of ~ 1.9 in ^'^^Nd, to ~ 3.3 
characteristic for the well deformed rotors ^^^Nd and ^^^Nd. For the = 90 isotope, in 
particular, the calculated i?4/2 = 3 and experimental R4/2 — 2.93 values, are very close to 
the characteristic, parameter-free, X(5) model prediction R4/2 = 2.91 [14]. Note, however, 
that our microscopic results for i?4/2, as well as other characteristic quantities, include the 
effect of /3 — 7 coupling, absent in the original X(5)-symmctry model. The panel on the right 
of Fig. 4 illustrates the evolution of another characteristic collective observable with neutron 
number: B(E2; 2^ — > Of) (in Weisskopf units). The important result here is not only that 
the calculation reproduces the swift increase of the empirical B(E2) values from less than 
30 Weisskopf units in ^^^Nd, to more than 160 Weisskopf units in ^^^Nd, but also that the 
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calculation is completely parameter-free. Namely, an important advantage of using structure 
models based on self-consistent mean-field single-particle solutions is the fact that physical 
observables, such as transition probabilities and spectroscopic quadrupole moments, are cal- 
culated in the full configuration space and there is no need for effective charges. Using the 
bare value of the proton charge in the electric quadrupole operator j(4{E2), the transition 
probabilities between eigenstates of the collective Hamiltonian can be directly compared 
with data. This is in contrast to, for instance, algebraic (interacting boson model) or shell- 
model approaches, that explicitly consider only valence nucleons and, therefore, calculate 
transition rates by adjusting the effective charges to reproduce some empirical values. For 
example, the X(5) predictions for transition strengths are normalized to the experimental 
B(E2; 2+ Of). We also note that the present resuffs, both for R4/2 and B(E2; 2+ Of), 
are in better agreement with data than those obtained in the axial GCM calculation with 
the Gogny force in Ref. [20]. In particular, the transition from spherical to well deformed 
Nd nuclei is less abrupt than in the axial GCM calculation and, therefore, the solutions of 
the collective Hamiltonian in five dimensions including /3 — 7 coupling, describe much better 
the transitional nuclei ^^^Nd and ^^°Nd. 

In Fig. 5 we display the isotopic dependence of the first excited 0^ state, and the ratio 
E{6~l)/E{02) in Nd nuclei. The microscopic excitation energies calculated with the PC-Fl 
energy density functional are compared with experimental values. The calculation repro- 
duces the empirical trend and, in particular, it predicts that the first excited O"*" state has 
the lowest excitation energy at = 90, in agreement with data. With the exception of 
the very low Oj state in ■'^^^Nd, the calculated excitation energies are also in quantitative 
agreement with experimental values. The fact that the band-head 0^ of the quasi-/? band 
has the lowest excitation energy in ^^°Nd can be attributed to the softness (fiatncss) of the 
potential with respect to P deformation (cf. Figs. 1 and 2). With the increase of the stiffness 
of the potential in ^^^Nd and ^^^Nd, the position of the /3-band is shifted to higher excitation 
energies, and the ratio E{6f)/E{02) decreases to less than 1/2. 

The prediction of a near degeneracy of the 6f level of the ground state band and the 
first excited 0+ state is another key signature of the X(5) model [5, 6]. In a very recent 
study of empirical order parameters for first-order nuclear QPT [36], Bonatsos et al. have 
shown that the ratio £^(6]'")/£^(0^) presents an effective order parameter of a first-order 
phase transition, and takes the special value of ~ 1 close to the phase transition point. 
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Furthermore, results of extensive IBM calculations with large boson number Nb show that 
the degeneracy of -E'(O^) and E{<o^) is a signature of a line of first-order transitions that 
extends across the Casten symmetry triangle. The ratios £^(6^)/i?(02 ) obtained with the 
PC-Fl energy density functional arc compared to the experimental values in Fig. 5(b). We 
notice that the calculation reproduces the abrupt decrease of this quantity between ^^^Nd 
and ^^^Nd, with a value close to 1 in ^^°Nd. The isotopic dependence of the theoretical 
values corresponds to the one predicted for a first-order phase transition in Ref . [36] . 

In Fig. 6 we compare the spectrum of the collective Hamiltonian for ^^°Nd with available 
data for positive parity states [34, 41, 42] , and with the predictions of the X(5) model. For the 
moments of inertia of the collective Hamiltonian we have multiplied the Inglis-Belyaev values 
Eq. (7) with a common factor determined in such a way that the calculated energy of the 2^ 
state coincides with the experimental value. In deformed nuclei, to a good approximation, 
the enhancement of the effective moment of inertia scales the relative excitation energies 
within each band by a common factor, but otherwise leaves the bandhead of the (5 band 
(0^) unaltered. Of course, this simple relation between the scaling of effective moments 
of inertia and excitation energies within each band would be exact only if rotational and 
vibrational degrees of freedom were decoupled. The degree of mixing between bands can 
be inferred from the distribution of iC-components (projection of the angular momentum 
on the body fixed-symmetry axis) of collective wave functions (cf. Table I). The transition 
rates are calculated in the full configuration space using bare charges. The inclusion of 
an additional scale parameter in our calculation was necessary because of the well known 
fact that the Inglis-Belyaev (IB) formula (7) predicts effective moments of inertia that are 
considerably smaller than empirical values. More realistic values are only obtained if one 
uses the Thouless-Valatin (TV) formula, but this procedure is computationally much more 
demanding, and it has not been implemented in the current version of the model. Here we 
rather follow the prescription of Ref. [37] where, by comparing the TV and IB moments of 
inertia as functions of the axial deformation for superdeformed bands in the A — 190 — 198 
mass region, it was shown that the Thouless-Valatin correction to the perturbative expression 
IB is almost independent of deformation, and docs not include significant new structures 
in the moments of inertia. It was thus suggested that the moments of inertia to be used 
in the collective Hamiltonian can be simply related to the IB values through the minimal 
prescription: Xk{q) — ll^{q){l + a), where q denotes the generic deformation parameter. 
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and q; is a constant that can be determined in a comparison with data. The value a ~ 0.4 
used for the excitation spectrum of ^^°Nd is comparable to those determined in the mass 
A = 190 - 198 region [37]. 

When the IB effective moment of inertia is renormalized to the empirical value, the exci- 
tation spectrum of the collective Hamiltonian determined by the PC-Fl density functional 
is in very good agreement with the available data for the ground-state band, (quasi) /? and 
7 bands. This is also true for the corresponding intraband and interband B(E2) values 
except, perhaps, for the strong transition 23 — > Oj 42 Weisskopf units), not seen in the 
experiment. Such a strong transition is probably due to the mixing between the theoretical 
(quasi) f3 and 7 bands. Tabic 1 includes the distributions of /^-components (projection of 
the angular momentum on the body fixed-symmetry axis), for the collective wave functions 
of the lowest three bands in ^^°Nd, ^^^Sm, and ^^^Gd. For ^™Nd, in particular, we notice 
a pronounced mixing of the states 2 J and 23 , assigned to the P and 7 bands, respectively. 
For higher angular momenta the mixing between the two bands is less pronounced. In ad- 
dition to the results for intraband transitions, we notice the very good agreement between 
the theoretical and experimental value B(E2; 0^ 2^). Again, we emphasize that the 
calculation of transition probabilities is parameter-free. Fig. 6 also includes the X(5) model 
predictions for ground-state {s — 1) and /?! (s — 2) bands of ^^°Nd. The theoretical spec- 
trum is normalized to the experimental energy of the state 2'1 and, in addition, the X(5) 
transition strengths are normahzed to the experimental B(E2; 2^ — > 0^). The simple X(5) 
model does not reproduce the data with the same accuracy as the solution of the collective 
Hamiltonian, especially the transition rates between bands. This could be due to the fact 
that -"^^^Nd is already slightly to the rotor side of the phase transition [38, 39], or simply to 
the fact that the X(5) model does not include — 7 coupling effects. Of course, even though 
the ansatz for the X(5) model involves a separation of variables in the P and 7 degrees of 
freedom, a full set of predictions for the quasi-7 band can be obtained [40]. However, this 
necessitates two additional parameters in the X(5) model that are adjusted, for instance, to 
the excitation energy of the quasi-7 band and the AK = 2 transitions. 

When compared with our particle-number projected axial GCM calculation of Ref. [22], 
where the same PC-Fl density functional was used, but configuration mixing included only 
prolate axially symmetric wave functions (cf. Fig. 2 of Ref. [22]), the full solution of the 
collective Hamiltonian in five dimensions displays reduced B(E2) values both for interband 
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and intraband transitions, and the (3 band is calculated at higher excitation energy. However, 
the inclusion of /3 — 7 couphng effects does not spoil the good agreement with data. The 
comparison with the X(5) model predictions is also illustrated in Fig. 7 where, for the yrast 
states of -"^"^^Nd, ^^'^Nd and -"^^^Nd, we compare the B(E2; L ^ L — 2) values and excitation 
energies calculated using the collective Hamiltonian based on the PC-Fl density functional, 
with the corresponding values predicted by the U(5) and SU(3) dynamical symmetries, and 
the X(5) model. Obviously the E2 rates and excitation energies for ^^°Nd are closest to those 
calculated from analytic expressions corresponding to the X(5) model. ^^^Nd does not differ 
very much from the X(5) limit, whereas the yrast states of ^^^Nd indicate that this nucleus 
is closer to a deformed rotor. We note, for instance, that in the case of the axial GCM 
model based on the Gogny interaction, the X(5) model predictions for the yrast excitation 
energies and transition rates were actually found to be closer to the calculated spectrum 
of ^'^^Nd, rather than ^^°Nd [20] . This comparison, of course, emphasizes the problem that 
the physical control parameter, i.e. the nucleon number, is not continuous and therefore in 
general a microscopic calculation cannot exactly reproduce the point of QPT, in contrast to 
phenomenological geometric or algebraic models that use continuous control parameters. 

A microscopic picture of the softness of the potential with respect to j3 deformation and 
the related phenomenon of QPT in ^^°Nd, emerges when considering the single-nucleon levels. 
In Figs. 8 and 9 we display the neutron single-particle levels in ^^^~^^^Nd, and the neutron 
and proton single-particle levels in ^^°Nd, respectively, as functions of the axial deformation 
parameter (5. The thick dot-dashed curves denote the position of the corresponding Fermi 
levels. In Fig. 8 we follow the evolution of the Fermi level with the increase in neutron 
number. For the range of deformation 0.2 < j3 < 0.4, in particular, the Fermi level crosses 
an energy interval of low level density. In ^^°Nd both the neutron and proton Fermi levels 
(cf. Fig. 9) are in the region of low level density for 0.2 < /3 < 0.4, and this result in 
the soft-/? potential (cf. Figs. 1 and 2). With further increase of the neutron number, the 
corresponding Fermi level in ^^^Nd approaches a region of high level density and the potential 
displays a pronounced, well deformed minimum at /? ~ 0.35. Further evidence for an abrupt 
change of structure between -"^^^Nd and -"^^^Nd can be clearly seen in the experimental spectra 
of the neighboring odd-Z nuclei: ^^^Pm and ^^^Pm [41]. Qualitatively very different band 
structures result from the couphng of the odd proton to the N — ^Q ox N — Q2 core nuclei. 
In ^^^Pm [N — 92) one finds regular A J = 1 low-lying rotational bands characteristic for the 
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limit of strong coupling. On the other hand, the lowest bands in ^^^Pm {N — 90) are regular 
AJ = 2 sequences, and this indicates that the odd proton is coupled to a core nucleus that 
is soft with respect to axial deformation. 

The microscopic results for the Sm and Gd isotopic chains are not very different from those 
of Nd nuclei, and thus in the remainder of this section we only present the most interesting 
features related to possible QPT. The self-consistent RMF+BCS triaxial quadrupole binding 
energy maps of I50.i52,i54gj^ -^^ ^j^g /3 — 7 plane (0 < 7 < 60^^) are shown in Fig. 10, and 
those of i52,i54,i56Q(^ pjg both cases we notice the increase of deformation of the 

prolate minima with increasing number of neutrons. However, in contrast to ^^°Nd, the two 
N = 90 isotones do not display pronounced /5-extended minima. This is especially the case 
for ^^^Sm, whereas a more extended fiat prolate minimum of the potential energy surface 
is calculated in ^^^Gd. Note also that the model does not predict the occurrence of oblate 
minima in the N — 92 isotones ^^^Sm and ^^^Gd. 

The slight change in the potentials of the two heavier = 90 isotones is also reflected in 
the corresponding spectra of the collective Hamiltonian. Figs. 12 and 13 display the spectra 
of ^^^Sm ^^^Gd, respectively, calculated with the PC-Fl relativistic density functional, in 
comparison with data [15, 35, 41, 42], and the X(5)-model predictions for the excitation 
energies, intraband and interband B(E2) values of the ground-state {s — 1) and /?i (s — 2) 
bands. For both nuclei the spectra of the collective Hamiltonian are in very good agreement 
with data, especially for the E2 transition rates. The agreement is perhaps not so good for 
the relative position of the /3 bands, which are calculated at somewhat higher excitation 
energies compared to data, and the effective moments of inertia of the quasi-7 bands are 
considerably smaller than the empirical values. As shown in Table 1, the mixing between the 
/3 and 7 bands is much less pronounced than in the case of ^^"^Nd. When compared to the 
predictions of the simple X(5) model, it is clear that the full collective Hamiltonian in five 
dimensions provides a more complete and accurate description of the low-energy spectra 
of ^^^Sm and ^^^Gd, especially for interband transitions between the P and ground-state 
bands. The comparison with X(5) is further illustrated in Figs. 14 and 15, where we plot 
the microscopic B(E2; L ^ L — 2) values and excitation energies of yrast states of Sm and 
Gd iV ?a 90 nuclei. While in the case of Nd isotopes the X(5) limit was very close to the 
results for ^^°Nd, for Sm and Gd the X(5) predictions lie in between the results obtained 
with the collective Hamiltonian for the = 88 and N — 90 isotopes. 
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IV. SUMMARY AND OUTLOOK 



The structural evolution of nuclei with neutron and/or proton number presents interesting 
shape transition phenomena, and in a number of systems signatures of first- and second- 
order quantum shape phase transitions have been observed. In this work we have applied 
the recently developed model for the solution of a five-dimensional collective Hamiltonian 
for quadrupolc vibrational and rotational degrees of freedom, with parameters determined 
by constrained self-consistent relativistic mean-field calculations for triaxial shapes, to a 
microscopic study of shape transitions in Nd, Sm and Gd isotopes with neutron number 
N Ki 90. Available data on excitation spectra and transition rates indicate that the A^" = 90 
isotones: ^^°Nd, ^^^Sm, and ^^^Gd, are located close to the point of first-order QPT between 
spherical to axially deformed shapes. 

Starting from self-consistent triaxial quadrupole binding energy maps in the /9 — 7 plane, 
calculated with the PC-Fl relativistic density functional, the parameters that determine 
the collective Hamiltonian: mass parameters, moments of inertia, and zero-point energy 
corrections, are calculated as functions of the deformations /3 and 7. The diagonalization of 
the resulting Hamiltonian yields the excitation energies and collective wave functions. An 
important feature of the model is that physical observables, such as E2 transition rates, 
are calculated in the full configuration space. Using the bare value of the proton charge, 
parameter-free transition probabilities between eigenstates of the collective Hamiltonian 
can be directly compared with experimental values. In the current implementation of the 
model the moments of inertia are calculated with the Inglis-Belyaev formula. The resulting 
values are considerably smaller than than the empirical moments of inertia and therefore, 
in order to compare the relative excitation energies with data, it was necessary to increase 
the Inglis-Belyaev moments of inertia. The effective values were adjusted to reproduce the 
experimental energies of the 2^ state in the ground-state band of each nucleus. 

A detailed analysis of shape transitions in Nd isotopes has been carried out, from the 
spherical ^^^Nd to the strongly deformed ^^^Nd. With increasing neutron number the Nd 
isotopes remain prolate deformed in the lowest state, and the calculation reproduces the 
characteristic signatures of shape transitions: R4/2 = E{A'^)/E{2'1) and B(E2; 2^ 0^), as 
well as the isotopic dependence of the effective order parameter E{6i)/E{02), which takes 
the special value of ~ 1 close to a first-order phase transition point. 
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Both the data and the results of microscopic calculations show that there is an abrupt 
change of structure ai N — 90. The microscopic potential energy surface of ^^°Nd displays 
a fiat prolate minimum that extends in the interval 0.2 < /3 < 0.4 of the axial deformation 
parameter, and a parabolic dependence on 7 for 7 < 30° in the region of the fiat prolate 
minimum. The resulting excitation spectrum of ^^*^Nd, calculated with the PC-Fl den- 
sity functional, and the corresponding intraband and interband B(E2) values, are in good 
agreement with available data for the ground-state band, (3 and 7 bands. The calculation 
reproduces the excitation energies of the bandheads 0^ and 2^ of the /3 and 7 bands, re- 
spectively and, when the Inglis-Belyaev effective moments of inertia are renormalized to the 
empirical value, also the relative excitation energies in all three bands. The parameter-free 
predictions for the B(E2)'s reproduce not only the experimental values for intraband transi- 
tions within the ground-state band and the band, but also for transitions from the /9 and 
7 bands to the ground-state band, including the measured branching ratios. The sequence 
of neutron single-particle levels in Nd isotopes, as functions of the axial deformation pa- 
rameter /3, and the isotopic dependence of the corresponding Fermi level, offer a qualitative 
explanation of the /3-softness of the potential in -^^^Nd. Similar results have been obtained 
for the Sm and Gd nuclei, and in particular the spectra of the collective Hamiltonian and 
the corresponding E2 transition rates reproduce the data in the A'" = 90 isotones ^^^Sm and 
is^Gd. 

A number of experimental studies over the last decade have disclosed candidate nuclei 
for quantum shape phase transitions in several mass regions. In addition to geometric and 
algebraic theoretical methods, it is also important to study these phenomena using micro- 
scopic models that explicitly take into account nucleonic degrees of freedom. Among the 
microscopic approaches to the nuclear many-body problem, the framework of nuclear energy 
density functionals (EDF) provides the most complete description of ground-state proper- 
ties and collective excitations over the whole nuclide chart. This work, together with similar 
recent studies, has shown that self-consistent mean-field models based on the EDF frame- 
work describe not only general features of shape transitions, but also particular properties 
of spectra and transition rates at the critical point of QPT. However, to calculate exci- 
tation spectra and transition probabilities, the self-consistent mean-field approach must be 
extended to include correlations related to restoration of broken symmetries and fluctuations 
of collective variables. This can be done either by performing GCM configuration mixing 



18 



calculations of projected wave functions, or by constructing collective Bohr-type Hamilto- 
nians with deformation-dependent parameters determined from self-consistent mean-field 
calculations. The possibility to perform self-consistent microscopic studies of shape transi- 
tions opens a new perspective on the origin of nuclear QPT in various mass regions. It is 
therefore important to systematically analyze, also employing different energy density func- 
tionals, various types of shape phase transitions that have been predicted in several regions 
of medium-heavy and heavy nuclei. 
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TABLE I: Distribution of the ii'-components K = 0,2,4 (projection of the angular momentum on 
the body fixed-symmetry axis) in percentage (%), for the conective wave functions of the lowest 
three bands in ^^^Nd, ^^^Sm, and ^^^Gd. 
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FIG. 1: (Color online) Self-consistent RMF+BCS triaxial quadrupole binding energy maps of the 
even-even ^^^"i^^Nd isotopes in the /3 — 7 plane (0 < 7 < 60*^). All energies are normalized with 
respect to the binding energy of the absolute minimum (red filled circle). The contours join points 
on the surface with the same energy (in MeV). 



23 




I I I I I I I I I I I I I I 

-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 

P 

FIG. 2: (Color online) Self-consistent RMF+BCS binding energy curves of the even-even ^^^"i^^Nd 
isotopes, as functions of the axial deformation parameter /?. Energies are normalized with respect 
to the binding energy of the absolute minimum for a given isotope. Negative values of /? correspond 
to the oblate (/? > 0, 7 = 60*^) axis on the /3 — 7 plane. 
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FIG. 3: (Color online) Self-consistent RMF+BCS binding energy curves of the ^^'^Nd nucleus, as 
functions of the deformation parameter 7, for two values of the axial deformation /3 = 0.2 and 0.25. 
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FIG. 4: (Color online) Evolution of the characteristic collective observables -R4/2 and B(E2; 2^ 
0^) (in Weisskopf units) with neutron number in Nd isotopes. The microscopic values calculated 
with the PC-Fl energy density functional are shown in comparison with data [34, 41, 42]. 
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FIG. 5: (Color online) Evolution of the first excited 0"*" state and the ratio E{6f)/E{02 ) with 
neutron number in Nd isotopes. The microscopic values calculated with the PC-FI energy density 
functional are compared to data [41, 42]. 
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FIG. 6: (Color online) The spectrum of ^^''Nd calculated with the PC-Fl relativistic density func- 
tional (left), compared with the data [34] (middle), and the X(5)-symmetry predictions (right) for 
the excitation energies, intraband and interband B(E2) values (in Weisskopf units) of the ground- 
state {s = 1) and Pi {s = 2) bands. The theoretical spectra are normalized to the experimental 
energy of the state 2f, and the X(5) transition strengths are normalized to the experimental B(E2; 
2+ ^ 0+). 28 




FIG. 7: (Color online) B(E2; L ^ L — 2) values (upper panel) and excitation energies (lower panel) 
for the yrast states in ^^^Nd, ^^'^Nd, and ^^^Nd, calculated with PC-Fl and compared with those 
predicted by the U(5), X(5), and SU(3) symmetries (open symbols). 
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FIG. 9: (Color online) Neutron (left) and proton (right) single-particle levels in ^^^Nd, as func- 
tions of the axial deformation parameter {3. Thick dot-dashed curves denote the position of the 
corresponding Fermi levels. 
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FIG. 10: (Color online) Self-consistent RMF+BCS triaxial quadrupole binding energy maps of 
I50,i52,i54gj^ in the /? - 7 plane (0 < 7 < 60°). All energies are normalized with respect to the 
binding energy of the absolute minimum (red filled circle). The contours join points on the surface 
with the same energy (in MeV). 
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FIG. 11: (Color online) Same as Fig. 11, but for the isotopes 152,154,156q^j_ 
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FIG. 12: (Color online) The spectrum of ^^■^Sm calculated with the PC-Fl relativistic density 

functional (left), compared with data (middle), and the X(5)-symmetry predictions (right) for the 

excitation energies, intraband and interband B(E2) values (in Weisskopf units) of the ground-state 

(s = 1) and Pi {s = 2) bands. The theoretical spectra are normalized to the experimental energy of 

the state 2f, and the X(5) transition strengths are normalized to the experimental B(E2; 2^ 0^). 

34 



1 ' 1 ' 1 — 7 1 ' — ~^ 




2 4 6 8 10 

Angular momentum (^) 



FIG. 14: (Color online) Same as Fig. 7, but for the isotopes I50=i52,i54gj^_ 
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FIG. 15: (Color online) Same as Fig. 7, but for the isotopes 152,154,156qj_ 
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